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The Muon {g — 2) experiment, E821, at the Brookhaven 
AGS has the goal to measure the muon anomalous mag- 
netic moment to a relative accuracy of ±3.5 x 10^^. A su- 
perferric 14 m diameter storage ring has been constructed 
and an averaged magnetic field uniformity over the 90 
mm diameter muon storage region of ± 1 part per million 
(ppm) has been achieved. A truncated double-cosine super- 
conducting septum magnet (the inflector) was constructed 
along with a fast non-ferric kicker. The performance of the 
storage ring, along with the physics results are reviewed. 

1 INTRODUCTION 

The measurement of magnetic moments has been im- 
portant in advancing our knowledge of sub-atomic physics 
since the famous 1921 paper of Stern, which laid out the 
principles of what we now call the "Stern-Gerlach experi- 
ment". The experimental and theoretical developments in 
the study of the electron's magnetic moment represent one 
of the great success stories of modern physics, with the ex- 
periment reaching a relative accuracy of ~ 4 parts in 10^ 
(parts per billion) [||] and the theory being constrained by 
our knowledge of the fine-structure constant a, rather than 
by the eight-order and tenth-order QED calculations. [Q| 

The gyromagnetic ratio g is defined by fls — g{e/2m)s, 
where s is the spin angular momentum, and /2 is the mag- 
netic moment resulting from this angular momentum. The 
Dirac equation predicts that g = 2, but radiative correc- 
tions increase the value at the part per mil level. The 
Particle Data Tables define the magnetic moment as /i = 
(1 + a)eTi/2m where a = — 2)/2 is the anomalous mag- 
netic moment (or simply the anomaly). 

When E821 began in the early 1980s, was known to 
7.3 parts per million (ppm).|^ The £82 1 Collaboration has 
reported three new measurements with relative accuracies 
of 13, 5 and 1.3 ppm respectively.!^ 0. ^ 

To the level of the experimental accuracy, the electron 
anomaly can be described by the QED of and photons, 
with the contribution of heavier virtual particles entering at 
a level below 4 ppb. The larger mass of the muon permits 
heavier virtual particles to contribute, and the enhancement 
factor is ~ {m^/mef ~ 40,000. The CERN measure- 
ment observed the effect on a^^ of virtual hadrons at the 
10 standard deviation level, [^l The standard model value 
of consists of QED, strong interaction and weak radia- 
tive corrections, and a significant deviation from the cal- 
culated standard model value would represent a signal for 
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2 THE EXPERIMENTAL TECHNIQUE 

The method used in the third CERN experiment and the 
BNL experiment are very similar, save the use of direct 
muon injection into the storage ring which was developed 
by the E821 collaboration. They are based on the fact that 
for g ^2 (or more precisely > 0) the spin gets ahead of 
the momentum vector when a muon travels transversely to 
a magnetic field. The Larmor and Thomas spin-precession 
and the momentum precession frequencies are 

geB eB eB 

ujs = '-x ^(1-7) — ; = — (1) 
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and the difference frequency gives the frequency with 
which the spin precesses relative to the momentum, 

uja = L^s - ^^c ^ i^;—) — (2) 
2 mc 

which is proportional to the anomaly, rather than to the 
full magnetic moment. A precision measurement of a^t re- 
quires precision measurements of the precession frequency 
uja and the magnetic field. The muon frequency can be 
measured as accurately as the counting statistics and detec- 
tor apparatus permit. The design goal for the NMR mag- 
netometer and calibration system was a field accuracy of 
about 0.1 ppm. The B which enters in Eq. ||is the average 
field seen by the ensemble of muons in the storage ring, 
< B >0=< J M{r,9)B{r,9)rdrd9 >^ where <j) is the 
azimuthal angle around the ring, r, 6 are the coordinates at 
a single slice of azimuth centered at the middle of the 90 
mm diameter muon storage region. M{r, 9) is the moment 
(multipole) distribution of the muon distribution, and cou- 
ples multipole by multipole with the magnetic field. It is 
very difficult to obtain adequate information on the higher 
moments of the muon distribution in the storage ring, so 
the presence of higher multipoles in the magnetic field is 
undesirable. 

The need for vertical focusing implies that a gradient 
field is needed, but the usual magnetic gradient used in 
storage rings is ruled out in our case. A sufficient mag- 
netic gradient for vertical focusing would spoil the ability 
to use NMR to measure the magnetic field to the necessary 
accuracy, and would also require detailed knowledge of the 
muon distribution. 

An electric quadrupole is used instead, taking advantage 
of the "magic" 7 = 29.3 at which an electric field does 
not contribute to the spin motion relative to the momen- 
tum. This can be understood from the famous Thomas- 
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which reduces to Eq. ^ in the absence of an electric field. 
Note that for muons with 7 = 29.3 in an electric field 
alone, the spin would follow the momentum vector 

The arrangement of a magnetic dipole field combined 
with an electric quadrupole field is called a Penning trap in 
atomic physics. However with a 14 m diameter and ~ 700 
T weight, the scale of our trap is quite different from the 
usual one. [|]] 

In order to meet the conditions discussed above, a goal 
of ±1 ppm uniformity of the < B > -field over the stor- 
age region was set and met. A round beam profile was 
chosen, since sharp corners would imply large higher mo- 
ments for M(r, 9). Given the projected knowledge of the 
muon distribution, the allowable strength of the quadrupole 
and higher magnetic multipoles was also determined. The 
quadrupoles are arranged in a four-fold symmetry cover- 
ing 43% of the ring. This geometry has the advantage that 
/?max — /3min to about 5% SO the beam does not breathe 
very much, making the average field calculation easier 

In our optimal running conditions, the 24 GeV/c proton 
beam in the AGS is accelerated in 12 proton bunches, each 
with an intensity of about 5 x 10^^ protons. The bunches 
are extracted one at a time at 33 ms intervals, and brought 
down a transport line to a Ni production target. The time 
distribution of the proton beam has at — 25 ns. The AGS 
cycle time is about 2.8 s. 

Pions at 0° are momentum analyzed and then brought 
into a 72 m straight FODO decay channel where muons are 
born. A second momentum slit permits one to choose for- 
ward muons about 1.6% below the pion momentum. These 
forward muons have > 90% longitudinal polarization, and 
are injected into the ring where their spin precesses. We 
store on the order of 10^ muons per fill of the ring. 

Unlike the conventional storage ring made of lumped 
elements, the requirement of a field uniform at the ppm 
level precluded breaks in the storage ring magnet. Thus 
the beam must be brought in through the fringe field of 
the storage ring to a point close the the central orbit. This 
is achieved by the use of a 1.7 m superconducting inflec- 
tor magnet[||| which nulls the field where the beam en- 
ters, but does not leak flux into the storage region, spoil- 
ing the precision field. An elevation view of the inflector 
exit, and the magnet is shown in Fig. |l]. This unconven- 
tional set of conditions means that there is no phase-space 
matching between the incoming beam and the storage ring. 
This mismatch reduces the calculated injection efficiency 
to - 8.7%. 

A kick of about 0. 1 Tm is needed to bring the beam onto 
a stable orbit. This is achieved with three 1.7 m long ferrite- 
free kickers, [p^ which can be thought of as single-loop 
pulsed magnets carrying a current of 4,200 A. The mini- 
mum inductance achievable of 1 .6 fiH limited the peak cur- 
rent to 4200A, and resulted in a current-pulse width ^2.5 
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Figure 1: The inflector exit-vacuum chamber geometry. 
The center of the storage ring is to the right. The gap be- 
tween the pole pieces is 180 mm, and the inflector exit is 
18 X 56 mm^ (ignoring the chamfer on the outer radius 
corners). 



times greater than the cyclotron period. This less than op- 
timal kicker pulse reduces the injection efficiency to about 
7.3%. Nevertheless the number of stored muons per fill is 
almost a factor of 100 over that available in the final CERN 
experiment,[|j and the injection-related background seen 
by the detectors is down by a factor of 50. [|[ 0] 

The precision magnetic field is at the heart of the 
experiment. [|ll| A profile of the storage ring magnet is 
given in Fig. 
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Figure 2: A schematic of the magnet cross section. 

The air-gap dominates the reluctance of the magnetic cir- 
cuit, and the pole-pieces are made of very low carbon, high 
quality magnet steel. The angle of the iron wedge in the air 
gap is used to eliminate the quadrupole inherent in a "C" 
magnet, and the wedge can be moved radially to adjust the 
dipole locally. The pole bumps were ground individually to 
minimize the sextupole component, and the pole face wind- 
ings (programmable current sheet) permit one to minimize 
the higher multipoles (on average), which do not vary much 
around the ring. 



Multipole expansion of B field in 1999 | 




-4 -3 -2 -1 A 1 2 3 4 
X (cm) 



Multipoles (ppm) 
normal skew 
Quad -2.23 2.17 
Sext -1.15 2.52 
Octu -1.33 1.85 
Decu 0.92 0.95 



Table 1: Parameters of the storage ring. 



Figure 3: Contour map of the magnetic field averaged over 
azimuth for 1999. The contours show 1 ppm differences. 



Parameter 


Value 


{g — 2) Frequency 
{g - 2) Period 


fa ^ 0.23 X 10*^ Hz 

Ta = 4.37/iS 


Muon Kinematics 


3.094 GeV/c 
7m = 29.3 
jT — 64.4 /is 


Cyclotron Period 


Tcyc = 149 ns 


Central Radius 


p = 7112 mm 


Magnetic Field 


Ba = 1.451 T 


Storage Aperture 


9.0 cm circle 


In one lifetime: 


432 revolutions 
14.7 (g-2) periods 



The success in shimming the magnet can be seen in 
Figs. 1^ and ^ which show the average field from 1999 and 
2000. The poorer field quality in '99 came from a dam- 
aged passive superconducting shield around the inflector, 
which permitted flux leakage into the storage ring. In 2000 
the damaged inflector was replaced, and we met the goal of 
±1 ppm uniformity. Several important parameters of the 
storage ring are given in Table |l]. 
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Figure 4: Contour map of the magnetic field averaged over 
azimuth for 2000. The contours show 1 ppm differences. 

The experimental signal is the from /i^ decay. Muon 
decay is a three-body decay, so the 3. 1 GeV muons produce 
a continuum of positrons (electrons) from the endpoint en- 
ergy down. Since the highest energy are correlated with 
the muon spin, if one counts high energy as a function 
of time, one gets an exponential from muon decay modu- 
lated by the [g — 2) precession. The expected form for the 
positron time spectrum is 



fit) = NfiC-'-'ll + AcOs{uJat + ^)] 



(4) 



However, a Fourier analysis of the residuals from this 
five parameter fit shows a number of frequency compo- 
nents which can only be understood after a discussion of 
the beam dynamics in the ring (see Fig. ph. 
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Figure 5: A Fourier Transform of the fit residuals from a 
5-parameter fit to the 1999 data set. 
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3 BEAM DYNAMICS 

2) ring is a weak focusing ring with 



n 



(5) 



where k is the electric quadrupole gradient. Several n - 
values were used for data acquisition: n = 0.137, 0.142 
and 0. 122. The horizontal and vertical betatron frequencies 
are approximately given by 



fx = fcVT 



0.929/c fy^ fcV^~ 0.37 fc 



(6) 

where fc is the cyclotron frequency and the numerical val- 
ues assume n ~ 0.137. 

The detector acceptance depends on the radial position 
of the muon when it decays, so any coherent radial beam 
motion will amplitude modulate the decay distribution. 
The principal frequency will be the "Coherent Betatron 
Frequency" 



/cBO = fc - fx = (1 - Vl - n)fc 



(7) 



which is the frequency a single fixed detector sees the beam 
moving coherently back and forth. An alternate way of 



thinking about this frequency is to view the ring as a spec- 
trometer where the inflector exit is imaged at each succes- 
sive betatron wavelength. In principle an inverted image 
appears at half a betatron wavelength, but the radial image 
is spoiled by the ±0.5% momentum dispersion of the ring. 
A given detector will see the beam move radially with the 
CBO frequency, which is also the frequency that the hor- 
izontal waist precesses around the ring. However, since 
there is no dispersion in the vertical dimension, the vertical 
waist is reformed every half a wavelength. The CBO fre- 
quency and its sidebands are clearly visible in the Fourier 
spectrum, and the vertical waist is just seen. A number of 
frequencies in the ring are tabulated in Table || 

Table 2: Frequencies in the (g — 2) storage ring for 



n ^ 0.137. 



Quantity 


Expression 


Frequency 


Period 


fa 


TT^ dllB 

?.Trm,r. P 


0.23 MHz 


4.37 IIS 


fc 


V 


6.7 MHz 


149 ns 


fx 


VI - nfc 


6.23 MHz 


160 ns 


fy 




2.48 MHz 


402 ns 


/cBO 


fc ^ fx 


0.477 MHz 


2.10 Ais 


/vw 


fc - 2fy 


1.74 MHz 


0.574 ^s 



The tune plane is shown in Fig.^ which shows resonance 
lines up to fifth order. Of the three n-values used for data 
collection, the n — 0.137 tune had a CBO frequency un- 
comfortably close to the second harmonic of {g — 2) (see 
Table ^ putting an AM sideband just by the {g — 2) fre- 
quency. This nearby sideband and has forced us to to work 
very hard to understand the CBO and how its related phe- 
nomena affect the value of uja from the fits to the data. The 
2001 data set was taken with the other two tunes, which 
substantially reduced our sensitivity to the CBO. 




Figure 6: The tune plane showing the three operating points 
used during our three year's of running. 



4 MONITORING THE BEAM PROFILE 

Three tools were available to us to monitor the muon dis- 
tribution. By studying the beam debunching after injection, 
one can gain information on the distribution of equilibrium 
radii in the storage ring. A wire chamber system located 
at one position around the ring permitted us to measure the 
trajectories of the and by tracing the trajectory back to 
the point where it is tangent to the ring, reproduces the de- 
cay position to within a few mm. The ring was equipped 
with two sets of scintillating fiber beam monitors which 
could be plunged into the storage region. Each set con- 
sisted of an X and y plane of seven 0.5 mm diameter scin- 
tillating fibers which covered the beam region. 

In Fig. ^ the signal from a single detector is shown at 
two different times following injection. The bunched beam 
is seen very clearly in the top figure, with the 149 ns cy- 
clotron period being obvious. The slow amplitude modula- 
tion comes from the {g — 2) precession. By 36 fis the beam 
has largely de-bunched. In Fig. [s] the inferred distribution 
of equilibrium radii is shown along with that obtained from 
a monte carlo tracking code. The agreement is seen to be 
good, and the measured distribution was used in determin- 
ing the average magnetic field and the radial electric field 
correction. 



e Time Spectrum: t = 6 |a,s 




e"*" Time Spectrum: t = 36 |-i s 




Figure 7: The time spectrum of a single detector soon after 
injection. The spikes are separated by the cyclotron period 
of 149 ns. 

The scintillating fiber monitors showed clearly the verti- 
cal and horizontal tunes as expected. Fig. ^jthe beam cen- 
troid motion is shown, both with the quadrupoles powered 
asymmetrically during scraping, and symmetrically after 
scraping. A Fourier transform of the latter signal shows 
clearly the expected frequencies. The traceback system 
also clearly sees the CBO motion. Additional details on 
beam dynamics will be available in Ref. [ |l2[ ] 




Figure 8: The distribution of equilibrium radii obtained 
from the beam de-bunching. 
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Figure 9: The beam centroid motion and its Fourier trans- 
form. 

A physics resuh with a total error of 1.3 ppm was ob- 
tained from our 1999 data set.[^ The result was originally 
2.6 standard deviations above the predicted standard model 
value, creating an outpouring of interest, and much spec- 
ulation that the first evidence for supersymmetry had been 
obtained. However, a sign mistake in a small piece of the 
theoretical hadronic contribution improved the agreement 
to 1.6 cr. We have an additional data set with /i+ which 
will have a total error of around 0.8 ppm, and also a /i^ 
data set which should also have a similar error, although 
the distribution between statistical and systematic errors is 
different for the two data sets. 

The theory value is also being refined further, and we 
hope that both theory and experiment will have presented 



new values by September 2002. 

5 CONCLUSIONS 

The muon {g — 2) collaboration has successfully built 
and operated a superferric storage ring with unprecedented 
field uniformity over such a large volume. All aspects of 
the apparatus, inflector, main magnet, kicker, electrostatic 
quadrupoles, detector system and electronics have achieved 
the design specifications. Interest in our new results re- 
mains high, and we look forward to improvements from 
our theoretical colleagues, which will improve the sensitiv- 
ity of the comparison to the standard model. 
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